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Abstract 

We study the motion of a single helical vortex in an unbounded, 
inviscid, incompressible fluid. The vortex is an infinite tube whose 
centerline is a helix and whose cross section is a circle of small radius 
(compared to the radius of curvature) where the vorticity is uniform 
and parallel to the centerline. Ever since Joukowsky (1912) deduced 
that this vortex translates and rotates steadily without change of form, 
numerous attempts have been made to compute these self-induced ve¬ 
locities. Here we use Hardin’s (1982) solution for the velocity field 
to find new expressions for the vortex’s linear and angular velocities. 
Our results, verified by numerically computing the Helmholtz integral 
and the Rosenhead-Moore approximation to the Biot-Savart law, are 
more accurate than previous results over the whole range of values of 
the vortex pitch and cross-section. We then use the new formulas to 
study the advection of passive particles near the vortex; we find that 
the vortex’s motion and capacity to transport fluid depend on its pitch 
and cross section as follows: a thin vortex of small pitch moves fast 
and carries a small amount of fluid; a thick vortex of small pitch moves 
at intermediate velocities and is a moderate carrier itself but it pushes 
fluid forward along the helix axis; and a vortex of large pitch, whether 
thin or thick, moves slowly and carries a large amount of fluid. 


1 Introduction 

The realisation that concentrated helical vortices are frequent features in 
natural and man-made flows has grown steadily since Parsons (1901) discov¬ 
ered that “a small spiral vortex existed just behind the tips of the propeller 
blades” used to generate thrust to propel a vessel. The theoretical study 
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Figure 1: A segment of a thin helical vortex. The vortex extends indefinitely 
in both directions and its centerline is a helix of pitch L and radius R lying 
on the surface of an imaginary supporting cylinder. 


of helical vortices has an even longer history. Here we highlight the earli¬ 


est results, but see Ricca (1994) for a detailed review of classical studies of 


the self-induced velocity of a single helical vortex in an inviscid, unbounded 
fluid (figure [^. [Kelvin (1880) found that one mode of vibration of a cylin¬ 
drical vortex is a wave that deforms the vortex axis into a helix of small 
radius and large pitch, and that this wave propagates along the axis with 
constant speed. Fitzgerald (1899) assumed that this result was valid for 


helices of any pitch and radius and speculated about the flow induced by 
the vortex: “there will be, on the whole, a flow along the inside of the spi¬ 
ral, but the motion of the fluid is complex.” Joukowsky (1912) showed that 


a helical vortex moves steadily without change of form and found that its 
self-induced velocity is approximately equal to the velocity of an osculating 
ring vortex; he further stated that any number of equal helical vortices sym¬ 
metrically arranged with respect to a common axis form a steadily moving 
arrangement. 

In the last two decades there has been a renewed interest on helical 


vortices (Ricca, 

1994; Mezic et ah, 1998 

Kuibin and Okulov 

, 1998; Boersma 

and Wood, 

1999[ Wood and Boersma, 

2001; Okulov, 2004) 

. Some of these 


works have concentrated on the binormal component of the self-induced 
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Figure 2: The self-induced velocity of a helical vortex, represented as the 
sum of either the azimuthal and axial components [uq and Uz^ respectively) 
or the tangential and binormal components {ut and respectively). 


Widnall, 1972 

) or the particle motion in the vicinity of the vortex 

Mezic et al. 

, 1998; Andersen and Brpns 

, 2014). 


velocity and have implicitly or explicitly dismissed the tangential component 
(see figure . This is unfortunate because accurate values of the full self- 
induced velocity are essential for a proper analysis of the vortex stability 
(see, e.g., 

(see, e.g.. 

The objectives of this paper are, first, to obtain expressions for the vor¬ 
tex’s linear and angular velocities {U and 12, respectively) that are accurate 
over the whole range of pitch values for vortices of small cross section (when 
compared to the radius of curvature of their centerline); and, second, to 
apply the new formulas in the study of particle motion in the flow induced 
by helical vortices. Section contains our theoretical calculations, which 
are based on the velocity field obtained by Hardin (1982); the verification of 
the new formulas for U and 12 by numerical computation of the Helmholtz 
integral and the Rosenhead-Moore approximation to the Biot-Savart law; 


and a comparison with previous results 

Joukowsky[ 

1912 

Da Rios 


1916 

Levy and Forsdyke, 1928 

; Widnall, 1972; 

Mezic et al., 1998; 

Okulov, 

|2004). 


In section we determine the vortex’s capacity to carry fluid by analysing 
the helical stream function of Hardin (1982) in a system that moves with 
the vortex. We present our conclusions in section 
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2 Self-induced motion of the vortex 


A helical vortex is a thin tube of infinite length whose centerline is a helix of 
uniform pitch lying on the surface of an imaginary supporting cylinder (see 
figure [^. The centerline of the vortex is given, in Cartesian coordinates, as 
follows: 


X — i?COS0, 

y = i?sin0, 
z = L9I2ti^ 

where 6 is the angle around the cylinder’s axis, R is the radius of the helix 
and L is the pitch of the helix. 

On the vortex’s circular cross-section the vorticity is uniform in magni¬ 
tude and direction, which is parallel to the centerline’s tangent. Note that 
both, the circular shape and the uniform vorticity, are leading-order approx¬ 
imations only: in a steady solution of the Euler equations the vorticity varies 
linearly with the distance to the centre of curvature and the cross-section is 
a small perturbation of the circular shape. The vortex thus has circulation 
r = Tra^cj, where a is the radius of the cross-section and lj is the magni¬ 
tude of the vorticity. The four parameters that uniquely define the vortex 
are then T, i?, L and a but, since the value of T merely changes the time 
scale of the evolution, we are left with two non-dimensional parameters: the 
vortex radius ol — a/R and the vortex pitch r = Lj^'KR. 


2.1 Analytical results 


Because of the geometry of the problem the self-induced velocity of the vor¬ 
tex can be expressed in either cylindrical components, pointing in 

the azimuthal and axial directions, respectively; or in natural components, 
pointing in the tangent and binormal directions, respectively (see fig¬ 
ure . We will use cylindrical coordinates instead of the now more popular 
natural coordinates because the translation speed of the vortex, ?7, equals 
the axial component of the self-induced velocity, whereas the angular 
velocity of the vortex, f], equals the azimuthal component divided by the 
helix radius, uq!R. 

Hardin ( 1982| ) found the velocity field produced by an infinitely-thin 


helical vortex. His solution, expressed in cylindrical coordinates, is divided 
in an interior field (valid inside the supporting cylinder) and an exterior 
field (valid outside the supporting cylinder). Here we reproduce only the 
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azimuthal and axial components needed for the calculation of the vortex 
motion: 

FT? 

— -Si{r, 4>) if r < i? 


U0{r,(p) = 


Uz{r,(t)) = < 
where (p = 0 — zjl, I — L/27r and 


irrl 

F FT? 


( 1 ) 




( 2 ) 


if r > i? 


Siir,(p) = cosm0 


m=l 

oo 


S2{r,(f>) = XI (^) 

m=l ^ ^ 

are Kapteyn series involving the modified Bessel functions and Im-> and 
their corresponding derivatives and 

Since we are dealing with vortices of finite cross section, it is possible to 
compute the vortex’s self-induced velocity by evaluating the velocity field at 
two diametrically-opposed points on the surface of the tubular vortex. The 
uniformity of the vorticity on the cross section guarantees that the average 
of these velocities is the actual velocity at the centerline, i.e. the velocity 
of the vortex. This approach has been previously used by Ricca (1994), 


Boersma and Wood (1999), and Okulov (2004) for the computation of the 


binormal component; here we will use it for the computation of the axial and 
azimuthal components. For simplicity, we choose to evaluate the velocity at 
points (r, 0 , 2 :) = (i?±a,0,0) or, in helical coordinates, (r, 0) = (i?±a,0). 
Therefore, the self-induced motion of the vortex is given by 


U = - a,0)-hi^z(i? + a,0)] 

n = ■^[u0{R-a,O) + ue{R + a,O)] 


( 3 ) 

( 4 ) 


A technical issue with this method is that the series Si and S 2 converge 
slowly, particularly when approaching the vortex itself, which hinders the 
computation of U and fl. This problem was solved by [Boersma and Wood 


5 


















Figure 3: The linear and angular velocities of a helical vortex [U and 
respectively) as functions of its pitch and radius (r and a, respectively). 
The thick contours indicate higher absolute values of ?7 or the grey areas 
indicate the region of the parameter space where is negative (clockwise 
rotation of the helix). 


(1999), who eliminated the singularities from Si and S 2 for the two particular 


points needed in the computation of the self induced motion: 


5i(i?-a,0) = (-^ + ln(e) + ln 

52(i? + a,0) = F + Me) + ln 





-WiT) + oil) 


where e = a/R{l + r^), o(l) represents an expression that tends to 0 when 
e ^ 0, and 


r 

W{t) = / 


sim t 


1 




_ + sin^ t)^/^ (1 + r^)^/^ 


dt 


with H{-) the unit step function. [Boersma and Wood (1999) computed 
W{t) by numerical quadrature for 21 values of r in the range [0.01 — 10] 
and obtained approximate analytical forms for large values of r (giving less 
than 1% error for r > 3) and for small values of r (giving less than 1% 
error for r < 0.4). Here we computed the value of W{r) using the numerical 


method described by Boersma and Wood (1999). 

The substitution of these forms of Si and S 2 in 0 -@ and the subse¬ 
quent substitution of the resulting expressions for uq and Uz in @-0 give, 
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after some algebra, the linear and angular velocities in nondimensional form: 

= (1 + ^2)3/2 - ln(VT + 7^) + (1 + T^f/^W{r)) 

^ /2(1 + r^) + ln(2/e) - ln(\/rT^) - (1 + - VK(t)) 

(1 + Y e^(l + — 1 


In dimensional form the self-induced motion of the vortex is given by 


U 

n 


AttR 

r 








( 5 ) 

( 6 ) 


Figure [ 3 ] shows U and f] in the region 10“^ < a < 0.4 and 0.1 < r < 10. 
All vortices translate in the direction of the axial component of the vorticity 
with a velocity U that increases as their pitch r and radius a decrease. The 
vortices generally rotate in a direction opposite to the azimuthal component 
of the vorticity, i.e. clockwise when seen from the direction in which the 
vortex translates. They rotate with an angular velocity Q that, for a fixed 
radius a, has a maximum absolute value when their pitch r is about one and 
that, for a fixed r, increases in magnitude as a decreases. Only relatively 
thick vortices of small pitch rotate in anti-clockwise sense (upper left corner 
in figure . 

Note that the translation speed U behaves as the helix’s curvature, 
l/i?(l + T^), which has a maximum value (1/i?) when r = 0 and then de¬ 
creases to zero as r goes to infinity. Similarly, the magnitude of the rotation 
speed behaves as the helix’s torsion, T/i?(l -h r^), which grows from zero 
when r = 0 to its maximum value (l/2i?) when r — 1 and then decreases 
to zero as r goes to infinity. 


2.2 Comparison with numerical results 

We verified equations (© and by computing the self-induced velocities 
of the vortex by numerical integration of the Helmholtz formula. 




X — X 


a; 


\x — X 


/|2 


dV, 


( 7 ) 


and the Rosenhead-Moore approximation to the Biot-Savart law. 


r /* [x — r(5)] X ds 

J {\x — 


( 8 ) 
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where a is the radius of the cross section of the tubular vortex and jjl is 
a parameter that depends on the vortex local structure, i.e. its vorticity 
distribution, curvature and torsion. Since, by our definition, a helical vortex 
has everywhere the same local structure /i is a constant but its value must 
be determined. This is usually done by choosing fi so that the integral 
(§ produces results in agr eement with a known solution obtained with a 
different method (see, e.g.. 


Saffman 


1995). Here we use fi — e because 


this value leads to the correct velocity for a ring vortex of uniform vorticity, 
constant curvature and zero torsion; this choice may thus be expected to 
give a good approximation for a helical vortex of uniform vorticity, constant 
curvature and constant torsion. 

In the computation of the Rosenhead-Moore integral we used a single 
filament represented by M nodes. In the computation of the Helmholtz 
integral we used 100 filaments, each represented by M nodes; none of these 
filaments coincides with the vortex centerline so, in contrast to the Biot- 
Savart case, there is no need to smooth out the singularity and, consequently, 
there is no free tuning parameter in Q. The value of M varied between 10^ 
and 10"^ depending on the helix pitch L. 

In order to compare the numerical and analytical results, we define the 
relative difference between them as follows 


Dj^ 


f {Un - U)^ + {mn - 

V U^ + 


1/2 


where quantities without subscript are obtained with formulas and 

those with subscript are obtained with either equation Q or ([8). We also 
find useful to define the relative difference between the linear velocities, 
Djj = \Un — U\/U ^ and the absolute difference between the angular velocities, 
Afl = (we must use this because is zero on a curve a — /(r), 

see figure and the relative difference would grow indefinitely even if the 
absolute difference were within round-off error). 

Figure shows Du^ Dq and Dt as functions of the vortex pitch (O.I < 
T < 10) for a given vortex radius {a = 0.05). In this case the velocities 
computed with (Ibj)-® and the Helmholtz integral have a mean relative 
difference Dt = 0.17% with a maximum of 0.27 % at r 0.65. In contrast, 
when the former are compared with Rosenhead-Moore calculations (using 
fi — e“^/^) the mean relative difference is 4.6% with a maximum of 5.7 % at 
T 1.4, the maximum Djj is 4.7% a,t r ^ and the maximum Afl occurs at 
T ^ 0.7. It is significant that the largest differences are obtained for vortices 
with high torsion, for it must be remembered that ji = is tuned to 

produce the correct velocity of a vortex of zero torsion. Therefore we varied 








Figure 4: The relative difference in linear velocity Djj^ the absolute differ¬ 
ence in angular velocity Af], and the relative difference in total velocity Dt-s 
as functions of the vortex pitch (r), for a given vortex radius {a = 0.05). 
Results obtained with equations [5][^ are compared with numerical integra¬ 
tion of equation]^ (continuous line), equation using /i = (dashed 

line), and equationusing fi = (dash-dotted line). 


the parameter ji in order to achieve a better agreement and found that 
ji — produces a mean relative difference Dt = 0.13% with a maximum 
of 0.25%. 

Identical results in Djj^ Dq and Dt are obtained if the Rosenhead-Moore 


kernel is replaced by the Winckelmans-Leonard kernel (Winckelmans and 
Leonard, 1993); i.e., if instead of equation Q we use 


^ ^ 2 ^ 2 ? 5/2 ^ 

J (|x — r{s)\^ -h iJ?o?) ' 

With this kernel it is /i = which leads to the correct ring-vortex 

velocity and 6% differences between analytically and numerically computed 
velocities of helical vortices, whereas ji = brings these differences down 
to 0.25%. 


2.3 Comparison with previous results 


J oukowsk^ ( 1912] ) an d |Da Rios| ( [TM^ found that, to the order of approxima¬ 
tion they used, the self-induced velocity is entirely in the binormal direction. 
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Figure 5: Linear and angular velocities, U and f], respectively, as functions 
of the vortex pitch (r), for a given vortex radius {a = 0.1). 


They found the binormal velocity, to be given by 

-—^ log ( —^ ) (Joukowsky 1912) 
AivRc \ a J 




^ logf- 


(9) 


AttR 


(Da Rios, 1916) 


where Rc = R{1 -\- is the radius of curvature. Note that the formula of 
Da Rios ( 1916[ ) is a straightforward application of the localised induction 
approximation he introduced ten years earlier (Da Rios, 1906). We decom¬ 


posed these velocities in cylindrical components in order to compute U and 
Q (see the thin blue lines in figure]^; they are in reasonably good agreement 
wit h ([^-([^, shown in black thick lines, for r > 2 only. 

Levy and Forsdyke ( 1928D and |Widnall (1972) computed the linear ve¬ 


locity U and the angular velocity Q of helical vortices of small pitch (in the 
ranges 0.25 < r < 1.25 and 0.1 < r < 1, respectively). [Levy and Forsdyfe 


(1928) avoided the singularity in the Biot-Savart law by evaluating part of 


the integral near, instead of on, the filament; Widnall (1972) obtained the 


self-induced velocity using the cut-off method and matched asymptotic ex¬ 
pansions. In the range 0.3 < r < 1 the results of Widnall (1972) and Levy 


and Forsdyke (1928), shown in dashed and dash-dotted lines in figure]^ are 


in better agreement with equations than any earlier or later result. 

Ricca (1994), Boersma and Wood (1999) and Mezic et al. (1998) argued 
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that the tangential component is unimportant and computed the binormal 
component only. It is true that the former contributes nothing to the time 
evolution of the Eulerian velocity field but it cannot be neglected in the 
analysis of particle motion (in section we will show that neglecting the 
tangential component, which can reach up to 30 % of the value of the binor¬ 
mal component for relatively thick vortices of small pitch, leads to significant 
errors in the determination of the flow topology and the capacity of the vor¬ 
tex to carry fluid). Mezic et al. (1998) decomposed the binormal velocity 


U}) to obtain U and the former shows reasonable agreement with (5) for 
all values of the vortex pitch, the latter shows good agreement with (6) for 
T > 1 only. 

The results of |Oknlov| (2004) for the binormal component not shown 
here) are in good agreement with ours; his results for U and fi, shown in 
green lines in figure differ greatly from ours as a consequence of his use 
of inconsistent assumptions, as explained below. In what follows the pre¬ 
fixes “O” and “R” indicate equations in lOkulovI (|2004|) and iRiccal (1 19941), 


respectively. Equation 04.4, which gives 14 as a function of r and was 
obtained using 02.4 and transformation relations between cylindrical and 
natural components (see R4.1). Then was computed under the assump¬ 
tion that the vorticity is uniform on the vortex cross section (see 04.6 and 
R4.12). A conflict exists because 02.4 implies 02.9, which means that the 
vorticity is parallel to helical lines of constant dimensional pitch (/ = L/ 27 r). 
Such a vorticity field cannot be uniform on the vortex cross section; thus 
04.4 gives an incorrect value of 14. The same conflicting assumptions are 


used by Okulov and Sprensen (2007, 2010). 


3 Flow topology 

With an appropriate change of variables it is possible to write the velocity 
field (ur^UQ^Uz) in terms of a stream function (Hardin, 1982). If we define 


(\) — d — zjl then U(^ — uq — ruzjl and the velocity field is given by 


Ur — - 


U(f) 


1 

r dcp 
dpj 
dr 
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with the stream function defined as follows (Hardin, 1982): 

r(r2 - i?2) rRr _ 


ip{r, (f)) = 


AttP 


ttP 


■S^ir^cf)) \ir<R 


r f ^Rf o / -f Tj 


where 


S3ir,(p) = cosm(j) 

m=l 
oo 

54(r,<^) = 


m=l 


fmr\ . / mR\ 


\ I) 


The stream function gives information about particle motion only if the flow 
is stationary. The flow induced by a helical vortex is time dependent in a 
reference frame fixed in space, but it can be made time independent in an 
infinite number of moving frames. The most natural choice, and the only 
one giving information about the capacity of the vortex to carry fluid, is a 
reference frame where the vortex is stationary; that is to say, in a frame that 
translates with linear velocity U and rotates with angular velocity fl. The 
steady stream function T is obtained with the simple transformation 

^ ^ i - y) (10) 


Because of the definition (/) = 0 — zjl^ the curves of constant T, which 
are streamlines of the helical flow represent the intersections of the 

stream surfaces of the three-dimensional flow {ur^ uq^Uz) with the polar plane 
z = 0. The intersections of the same stream surfaces with the meridional 
planes 0 = 0 , tt can be readily obtained using the same definition. 

The flow topology, or phase portrait in the language of dynamical- 
systems, consist of the set of stagnation points plus the streamlines that 
divide the flow in regions of qualitatively different behaviour. In the co¬ 
moving frame the vortex, located at (r, 9) — (a, 0), corresponds to a stagna¬ 
tion point of elliptic type and the symmetries of T indicate that other fixed 
points, when they exist, are located on the line 0 = 0,7r. We found these 
stagnation points using a numerical method in order to perform a systematic 
exploration of the region of the parameter space defined by 10' < a < 0.4 

and 0.1 < r < 10. 
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Figure 6: The three flow regimes in the parameter plane (r, cr): (I) large- 
pitch helices; (II) thin small-pitch helices; and (III) thick small-pitch helices. 
The thick lines show the boundaries found in this work, the thin lines show 
those found by Mezic et al. (1998) and the dashed lines show the results of 


Andersen and Brpns (2014). 
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Figure 7: The helical stream function T on the polar plane (r, 9) for repre¬ 
sentative cases of each regime: (a) Large-pitch helices (t=0.8, (a=0.05) , (b) 
thin small-pitch helices (t=0.3, (a=0.0001), and (c) thick small-pitch helices 
(t=0.2, a^Om). 


We found three qualitatively different ffow topologies (figure [^. These 
were first identified by Mezic et al. (1998) and Andersen and Brpns (2014); 
they, however, only took into account the binormal component of the vortex 
Mezic et al. ( 1998[) c omputed using the results of Ricca| ( |1994 ) 
Andersen and Br0ns| ( |2014| ) used equation ^ with log{Rc/a). As 


motion: 

whereas 


shown by the blue lines in figure all calculations that do not include the 
tangential component give large errors in U and fl, particularly for small 
values of the vortex pitch (r < 1). Consequently Mezic et al.| ( |1998 ) and 
Andersen and Brpns ( 2014[ ) found regime boundaries that are significantly 


shifted from the ones found here (see the thin lines in figure |^. 

The properties of each ffow topology, and the consequences for the vor¬ 
tex’s capacity to carry fluid, are the following: 

Regime I (which occupies mostly r > 1 for all values of cr, see figure]^ 
is characterised by the existence of only one stagnation point in the plane 
r—0, this is of elliptic type and coincides with the vortex itself (see figure]^). 
The presence of only one elliptic stagnation point implies that all particles 
have, qualitatively, the same behaviour: they all rotate around the vortex. 
There is, however, an important quantitative difference: while the particles 
that are close to the vortex drift slowly downwards, the particles that are 
farther from it drift with a speed that rapidly approaches —U (figure]^). 
In the fixed frame this means that nearby particles follow the vortex in its 
translation while distant particles remain close to the horizontal plane where 
they were initially located. 

Regime II (which occupies roughly r < 1 and a < 0.003, see figure]^ is 
characterised by the existence of two elliptic stagnation points and a hyper¬ 
bolic stagnation point between them (see figure]^). Now there are three 
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Figure 8: Same as figure but now on the meridional planes (r, z) corre¬ 
sponding to 0 = 0, TT. 


regions with qualitatively different behaviour: one set of particles rotate 
around the vortex and are permanently trapped by it; a second set of parti¬ 
cles rotate around a line close to the helix axis while drifting downwards; and 
the exterior particles rotate around those two sets while drifting downwards 
(see figure [^). 

Regime III (which occupies roughly r < 0.4 and a > 0.0001, see fig¬ 
ure]^ is characterised by the existence of two elliptic stagnation points and 
a hyperbolic stagnation point next to them (see figure]^). There are again 
three regions with qualitatively different behaviour: one set of particles ro¬ 
tate anti-clockwise around the vortex and are permanently trapped by it; 
a second set of particles rotate clockwise around a line close to the helix 
axis while drifting upwards; and the exterior particles rotate anti-clockwise 
around those two sets while drifting downwards (see figure |^). This is the 
type of flow Fitzgerald (1899) speculated about. 
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4 Conclusions 


We have found new expressions for the linear and angular velocities of helical 
vortices that have, to leading order, uniform vorticity and circular cross- 
section. These expressions are valid for vortices of any pitch (r) when their 
cross-sectional radius (a) is much smaller than their radius of curvature (i.e., 
cr <C 1 + T^). 

Numerical computations of the motion of helical vortices using the Helmholtz 
integral and the Rosenhead-Moore approximation to the Biot-Savart law 
show that the new formulas give the correct velocities or, at least, very good 
approximations. Indeed, if we assume that ji — in the Rosenhead- 

Moore kernel leads to the correct velocity then the best approximations are 
given by the new formulas, within 6 


Widnall (1972), 13.6 %; and Levy 


and Forsdyke (1928), 13.9 %. If, on the other hand, we assume that it is 


fi — e~^ that leads to the correct velocity then the best approximations 
are given by the new formulas, within 0.5 %; Widnall (1972), 14.1 %; and 
Mezic et al. (1998), 17.8 %. In any case, equations (IT® give the best 


approximation published so far. 

The pitch and radius of the vortex determine its motion and capacity 
to carry fluid in the following manner: large-pitch vortices, whether thin 
or thick, translate slowly and carry with them a large volume of fluid; thin 
small-pitch vortices translate fast and take with them a small volume of fluid; 
and thick small-pitch vortices translate at intermediate velocities, take with 
them a moderate volume of fluid and, more importantly, push fluid forward 
along the helix axis. 
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